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ABSTRACT 

The inversion of gravitational lens systems is hindered by the fact that multiple mass 
distributions are often equally compatible with the observed properties of the images. 
Besides using clear examples to illustrate the effect of the so-called monopole and mass 
sheet degeneracies, this article introduces the most general form of said mass sheet 
degeneracy. While the well known version of this degeneracy rescales a single source 
plane, this generalization allows any number of sources to be rescaled. Furthermore, 
it shows how it is possible to rescale each of those sources with a different scale factor. 
Apart from illustrating that the mass sheet degeneracy is not broken by the presence 
of multiple sources at different redshifts, it will become apparent that the newly con- 
structed mass distribution necessarily alters the existing mass density precisely at the 
locations of the images in the lens system, and that this change in mass density is 
linked to the factors with which the sources were rescaled. Combined with the fact 
that the monopole degeneracy introduces a large amount of uncertainty about the 
density in between the images, this means that both degeneracies are in fact closely 
related to substructure in the mass distribution. An example simulated lensing situ- 
ation based on an elliptical version of a Navarro-Frenk- White profile explicitly shows 
that such degeneracies are not easily broken by observational constraints, even when 
multiple sources are present. Instead, the fact that each lens inversion method makes 
certain assumptions, implicit or explicit, about the smoothness of the mass distribu- 
tion means that in practice the degeneracies are broken in an artificial manner rather 
than by observed properties of the lens system. 
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1 INTRODUCTION 

Not only can strong gravitational lensing yield impressive 
images, from multiply imaged quasars over deformed galax- 
ies to partial or full Einstein rings, it is also an invaluable 
tool to estimate the mass of the deflector. Moreover, since 
the precise positions and deformations of the images of a 
source depend on the exact distribution of the mass, gravi- 
tational lensing even holds the promise of constraining the 
shape of said mass distribution. 

In principle, this strong lens inversion sounds fairly 
straightforward: use a particular model and optimize its 
parameters so that it reproduces the observations as well 
as possible. In practice however, one is hindered by grav- 
itational lensing degeneracies which allow a wide range of 
mass distributions to produce the exact same image con- 
figuration. On the level of parametric lens inversion re- 
sults, this can manifest itself as p arameter degeneracies (e.g. 
IWambsganss fc PaczvnskH |l994|)), but at their core the de- 
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generacies are or course present at the level of t he mass dis- 
tribut ion . There, the mass sheet d egeneracy (e.g. lFalco et alj 
|l985t) or lSaha fc Williams! (120061 ) ) will rescale the projected 
mass distribution while adding a constant density mass sheet 
(or disc), only affecting the time delays between the im- 
ages. Recently it was shown that a similar procedure is 
still possi ble when multiple sources at different redshifts are 
observed (|Liesenborgs et al1l2008ah . While this degeneracy 
therefore can still be broken when time delay measurements 
are present, the situation is f ar worse with respect to t he so- 
called monopole degeneracy |Liesenborgs et al. 2008b). This 
type of degeneracy can be used to redistribute the mass in 
between the images, without changing any of the observable 
properties of the images. 

By presenting a further generalization of this mass sheet 
degeneracy, this article will illustrate an important relation 
between different degenerate versions of a lensing mass dis- 
tribution. This will illustrate how degeneracies are closely 
related to the substructure in the mass density, which, in 
turn, offers insights into what one may constrain about the 
lensing mass. 
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After briefly introducing the lensing formalism and used 
notations in section [2] the basics of the monopole and mass 
sheet degeneracies are reviewed in section [3] This will serve 
as the starting point for some extensions in section [4] after 
which an experiment in section[5]will illustrate precisely how 
difficult it can be to estimate the parameters of a specific 
type of mass distribution, an elliptical generaliz ation of a 
Navarro- Frenk- White profile (iNavarro et al.l I19961 ) . Finally, 
in section [U] we shall discuss what the implications of these 
degeneracies are and what different types of information can 
help constrain. 



2 LENSING FORMALISM 



In this section, the same notation as ISchneider etafl (|l992h 
is used, which the interested reader can consult for a thor- 
ough review. The lens equation ([T} in essence states that 
due to the gravitational deflection d of a light ray, when 
looking in a direction one sees what would be observed in 
direction f3 if the lens effect could be turned off: 



(i) 



The deflection angle a depends on the two-dimensional pro- 
jected mass distribution E in the following fashion: 
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(2) 



In these equations, the geometry of the lensing scenario is 
given by the angular diameter distances D d , D s and D ds , 
describing the distance to the lens, to the source and the 
distance between lens and source respectively. It will often 
be convenient to use the so-called critical density: 



E C r 



D. 
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(3) 



which depends on the redshift of a source through the an- 
gular diameter distances involved. 

In the special case of a circularly symmetric mass dis- 
tribution, for convenience assumed to be centered on the 
origin of the coordinate system, the expression of the deflec- 
tion angle reduces to: 

in which M(6) represents the mass enclosed within radius 0: 
e 

2 



M{ff) = 2-kDI I £(0')0'd0' 
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The deflection angle can be shown to be proportional 
to the gradient of the so-called lens potential 



(6) 



which is related to the time light takes to travel from the 
source at position /3 to an image at position 0: 

t{d)= l + z d D,D s ^ (0 _ pf ^ + congt (7) 

In this last equation, z d is the redshift of the gravitational 
lens. Since, in practice one is only interested in the differ- 
ence in time delay between images of the same source, the 



constant in this expression is of no importance. The curva- 
ture of the lens potential can be shown to represent the mass 
density at the corresponding point: 



(8) 



Finally, the magnification ^ that an image at position 
experiences, can be calculated from the lens equation as 
follows: 



dpi 



m _1 (0) 



i.e. it is the jacobian determinant of this equation. 



3 DEGENERACY BASICS 
3.1 Monopole degeneracy 



(9) 



The s o-called monopole degeneracy, first mentioned in Saha 



(2000) and studied in more detail in iLiesenborgs et al 



(2008b), is both the least complex and most elusive degener- 
acy one can think of. Consider the mass distribution shown 
in the inset in Fig. [T] Being a circularly symmetric mass dis- 
tribution, its lens effect is governed by equation Q. Since 
this particular mass distribution was constructed to have 
zero enclosed mass outside a particular radius, the same 
equation shows that it will not cause a deflection angle out- 
side said radius. This kind of mass density will be referred 
to as a monopole, since it can be described by only the 
monopole term in a multipole expansion. 

An example will illustrate the importance of such a 
monopole. The elliptical mass distribution shown in the left 
panel of Fig. [T] transforms a single source into five distinct 
images, as can be seen in the left panel of Fig. [2] Indicated in 
the same figure, is a circular region in which none of the im- 
ages are located, and which has the precise same diameter as 
the aforementioned monopole. When adding the monopole 
in this location, this implies that the lens equation will not 
be affected outside the circular region; only within will the 
deflection of light be affected. That this is indeed the case, 
can clearly be seen in the right panel of Fig. There, the 
lens effect is calculated for the same source, but for the com- 
bined mass distribution shown in the right panel of Fig. [T] 
As was expected, the exact same five images are generated 
by this clearly different mass density. The fact that the lens 
equation was modified inside the circular region can be seen 
in the change in critical lines in this region. Since the deflec- 
tion is different inside the circular area, in principle it is even 
possible that extra images are generated there. However, this 
particular example shows that this does not necessarily need 
to be the case, even for a rather large modification of the 
mass distribution. 

This example already shows that by 'borrowing' some 
mass from a relatively large region, it is easy to introduce 
(or remove) a mass peak in a mass density, without changing 
the images that are gene rated. As was explicitly shown in 
l|Liesenborgs et~aL 2008b) however, it is also possible to use 
such monopoles as basis functions which allows a complex 
redistribution of the mass by adjusting their weights. As 
long as no basis function overlaps with an image, none of 
the images will be affected by this operation. Some care 
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Figure 1. Left panel: this non-singular isothermal ellipse based gravitational lens transforms a single source into the five images that 
can be seen in the left panel of Fig. [2] Right panel: when the mass distribution shown in the inset is added in the region enclosed by a 
circle in Fig. [2] (left panel), this mass distribution is created. Note that although the added mass distribution itself contains a region of 
negative density, the total density of the lens is positive everywhere. 
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Figure 2. Left panel: the gravitational lens shown in the left panel of Fig. [T]creates five images of the source enclosed by a small square. 
A circle encloses a region which does not contain any of the images, and can therefore be used to illustrate the monopole degeneracy (see 
text). Right panel: when the circularly symmetric mass distribution from the inset of Fig. [T] is added in the circular region, the resulting 
mass distribution (right panel of Fig. [TJ still maps the source onto the exact same images. Note that the change in the critical lines and 
caustics illustrates the fact that inside the circular region, the lens equation indeed has changed. 



needs to be taken not to introduce additional, unobserved 
images though. 

From the discussion above, it is already trivially clear 
that this type of degeneracy does not modify the image po- 
sitions generated by a particular source. The fact that the 
deflection angle is zero outside a certain radius, by equa- 
tion © implies that the magnification at the location of 
the images will not be affected either. Equation (|6} on the 
other hand shows that the lens potential outside said radius 
may change by a constant value. However, when measuring 
time delays between images, such a constant offset in the 
potential will have no effect, as can be seen in equation ([7]). 



3.2 The mass sheet degeneracy 

The mass-sheet degeneracy is undoubtedly one of the most 
famous degeneracies in gravitational lensing. Contrary to 
what the name may suggest, to obtain a degenerate solu- 
tion it is not sufficient to simply add a sheet of mass to an 
existing mass distribution. Instead, one has to rescale the 
original mass distribution as well, prompt ing the use of the 
altern ative name of steepness degeneracy l|Saha fc Williams! 
l200rj ). 

As a concrete example, let us call the mass distribution 
of the left panel of Fig.[T]So, producing a deflection field do. 
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Figure 3. The mass sheet degeneracy in its simplest form causes 
a rescaled version of the source plane to correspond to the same 
images. The fact that the entire source plane has been rescaled 
is most obvious when comparing the caustics to those of the left 
panel of Fig. [2] In this particular example, a value of A = 0.7 
was used in the construction of a degenerate solution according 
to equation (JTTJ . 



Using equation one can show that for a sheet of mass of 
precisely the critical density E cr , the deflection angle simply 
becomes: 



(10) 



If one then constructs a new mass density Ei as follows: 

S 1 (6>) = AE (6») + (l-A)E cr , (11) 

it is a straightforward exercise to show that for the combined 
deflection Adio + (1 — X)a a the new lens equation is: 



j3 1 (6) = A/3 (0) , 



(12) 



i.e. merely a rescaling of the original one. This means that a 
rescaled version of the source plane corresponds to the same 
images, as is illustrated in Fig. 

If the source involved is variable, it may be possible to 
measure time delays between the images: 



Atij=t(O i )-t(0 j ), 



(13) 



where Oi and 8j are the positions of two of the images. For 
a sheet of mass consisting of precisely the critical density, 
the projected potential is 



MO) = ~e 2 



(14) 



so that the mass sheet degeneracy transforms an initial lens 
potential ipo into 



Vi(0) = At/>o (<?) + ■ 



A. 



(15) 



Noting that the mass sheet degeneracy changes a source 
position (3 to /3 1 = A/3 , a brief calculation shows that the 
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Figure 4. In this example, two sources at a different redshift 
each produce two images due to the presence of a gravitational 
lens. Because of the different redshifts, each source would require 
a different critical density in the construction of a solution that 
differs by the mass sheet degeneracy. The enclosed mass profiles 
corresponding to these two densities are shown as curves Ma an d 
Mb- The locations of the images of each source are indicated on 
these profiles. The dashed line, an interpolation between the two 
curves, shows a total mass profile that effectively behaves as the 
two sheets, from the perspective of the images (see text). 



relationship between the original time delay Atij,o and the 
time delay of the degenerate version Atij,i becomes: 



At 



AAt 



(16) 



This means that time delay measurements break this mass 
sheet degeneracy, since a particular version of the degener- 
acy corresponds to a particular time delay. 

The mass sheet degeneracy also has a particularly sim- 
ple effect on the magnification factor /x. Since each dimension 
is scaled by a factor A, the surface area of a small source is 
scaled by a factor A 2 . Keeping the sizes of the images con- 
stant, this means that the new magnifications of the source 
are given by /ii = \~ 2 fio. This result can also be obtained 
by combining the general expression for the magnification 
(O and the relation between the degenerate versions of the 
lens equation (I12[) . 



4 GENERALIZED MASS SHEET 
DEGENERACY 

A first gen eralization of the mass sh eet degeneracy was de- 
scribed in iLiesenborgs et all i|2008al ). A brief review shall 
be presented here, before discussing an even more general 
version of this degeneracy. It has often been claimed that 
having two distinct sources at different redshifts breaks the 
mass sheet degeneracy. From the construction above, the 
reason for this is clear. A sheet with precisely the critical 
density was used, but since the critical density is redshift 
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dependent, each source would now require the use of a sheet 
with a different density. 

This situation is illustrated in Fig. [4] Source A would 
require a sheet of mass which gives rise to the enclosed mass 
profile Ma(9)- Similarly, source B would need a sheet with 
mass profile Mb{0). The radii at which images of each source 
are present are indicated on the corresponding mass profile. 
While it is true that a single mass sheet can no longer be used 
for a construction like the one in equation {TTJ, it indeed is 
possible to construct a mass density that for each source has 
the same effect as its corresponding mass sheet. 

The key point is that for a circularly symmetric mass 
distribution, equation Q shows that the deflection angle is 
caused by the total enclosed mass within a specific radius, 
without any dependence on the structure of the mass profile 
within said radius. This implies that if one were to construct 
the mass profile indicated by the dashed line in Fig. [4l the 
deflection angle at the location of the images of source A will 
be the same as when caused by a mass sheet of precisely 
its critical density. Similarly, the images of source B will 
effectively 'see' a mass sheet with their critical density. The 
mass density corresponding to this new profile can then be 
used to construct a degenerate solution in the same way as in 
equation (| llf) . Calling Ra and Rb the radii at which images 
of sources A and B are located, by construction this implies 
that: 

Ei(-Ra) = AS (i?A) + (l-A)E cr (2 A ) 

Ei(i? B ) = AE (i? B ) + (l-A)E cr (z B ). (17) 

The effect is the same as for the original mass sheet 
degeneracy: the images will now correspond to sources which 
are rescaled by the same factor A. The magnifications are 
therefore again modified by A~ 2 . The effect on the time delay 
is not as straightforward however. For a circularly symmetric 
lens, combining equations ([6]) and ((4]) shows that 

m k C ^ d "' (i8) 

i.e. the lens potential, and therefore the time delays, will 
depend on the precise manner in which the interpolation 
between the two curves in Fig. [4] is done. 

Using the generalization of the mass sheet degeneracy 
that was just shown, the sources are still rescaled with the 
same factor. But it is also possible to construct a degen- 
erate solution which allows one to scale each source with 
a different factor. To illustrate this, consider the situation 
shown in the left panel Fig. [5] an elliptical mass distribution 
causes two sources to produce three images each. In what 
follows, only source A will be rescaled with a certain factor. 
In a subsequent step the other source can be rescaled in a 
similar way, using a different factor. 

Focusing on only one source and temporarily forgetting 
about the other, the source can be rescaled using the classic 
mass sheet degeneracy from equation {TTJ. This means that 
to the initial Eo, a mass distribution E a dd was added to 
obtain Ei: 

E add (0) = (1 - A) [E cr (z A ) - E o (0)] . (19) 

Taking the second source back into account, it is clear that 
adding E a dd to the original mass distribution will cause the 
images of source A to correspond to a scaled version of that 
source, but the effect on the second source is less desirable. 



However, suppose it is possible to modify E add so that 
the deflection it causes for source A is unaltered, and no 
deflection is caused at the location of the images of source B. 
In that case, after adding E ad d, the images of source A would 
correspond to a scaled version of that source, and the images 
of source B would still be projected onto the same source 
area. Calling I\ and 7b the points of which respectively 
the images of sources A and B consist, a add (/B) = by 
equations © and ([5| implies that E ad d(-?B) = 0. If one then 
constructs Ei = Eo + E add , in general one finds that 

Ei(/a) = AEo(Ja) + (1 - A)E cr (^ A ) 

Ei(/ B ) = E (/b) (20) 

Starting from Ei, in a next step one can then exchange the 
roles of the two sources and rescale source B with a different 
factor. 

There still remains the question of how to modify E a dd 
so that it does not produce a deflection at the location of 
the images of source B. This can be done by using a large 
number of the monopole basis functions we met earlier. If 
none of these overlap with the images of source A, the de- 
flection there will be unaltered. To make sure that the de- 
flection angle at the location of the images of source B van- 
ishes, appropriate weights of these basis functions must be 
sought. An example of this approach can be seen in the 
center panel of Fig. [5] In this example, the weights were de- 
termined by a genetic algorithm, an optimization strategy 
inspi red by natural sel e ction. A similar procedure is used 
as m iLiesenborgs et al l l |2008bh : a trial solution consists of 
a set of weights for the monopole basis functions. Starting 
from randomly initialized trial solutions, the algorithm then 
'breeds' increasingly better trial solutions by combining ex- 
isting solutions and introducing mutations. A key step in 
this procedure is the application of selection pressure: bet- 
ter solutions should produce more offspring. In this case, a 
trial solution is deemed better than another simply when its 
deflection field differs less from the desired one at the loca- 
tions of the images of source B. Since, when generating a 
degenerate solution in this fashion, one still uses a rescaled 
version of the original mass distribution, this procedure still 
produces a mass map which appears rescaled, as can be seen 
in the center panel of Fig. [S] 

Another approach which does not cause such an ap- 
parent rescaling, can be used as well. From the discussion 
above, it is clear that the final version of the mass distribu- 
tion E add should produce a specific value of the deflection 
angle a add at the location of the images of source A. On 
the other hand, at the location of the images of source B, 
no deflection should be produced. We therefore know the 
deflection angle at specific locations, and at other locations 
the deflection angle can take a large number of values, the 
only real constraint being that no additional images should 
be produced. Interpolating a deflection field in such a way 
that V x a a dd = 0, or equivalently that the deflection orig- 
inates from a scalar lensing potential, is precise ly what is 
done in the LensPerfect lens inversion procedure l|Coe et al.l 
2008). Using a similar interpolation procedure to determine 
£*add and therefore E a d d , yields the result shown in the right 
panel of Fig. [5] While more difficult to ensure that the over- 
all mass distribution is positive, in this example it produces 
an interesting result. Comparing the resulting mass distri- 
bution with the original one in the left panel, one sees that 
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Figure 5. Left column: this gravitational lensing situation is used as a starting point to illustrate the general mass sheet degeneracy in 
which different sources are scaled with different values. The top figure depicts the elliptical mass distribution of the gravitational lens; 
the two other figures show how this lens creates three images of two different sources. In this example, only the source in the middle panel 
(source A) will be rescaled, the source in the bottom figure (source B) will not. Center column: after modifying the mass distribution in 
the left column using the method of the monopole basis functions (see text), a rescaled version of source A is needed to create the same 
images as before. Source B on the other hand does not need to be altered to recreate the same image configuration. To illustrate that 
the source plane of source A has been rescaled toward the origin can be seen by comparing the new source position to the original one. 
Right column: similar to the center c olumn, but in t his case a similar interpolation method for the deflection field was used as in the 
LensPerfect lens inversion procedure llCoe et al.ll2008Tl . 



the differences are rather subtle. And since no rescaling of 
the mass density itself seems to be involved, one would not 
suspect a mass sheet-like degeneracy to be at work here. 



Note that in principle, a similar construction is possible 
when only using point sources. When the deflection angles 
are only conserved at very specific point locations and not 
over a more extended area, the gradients of the deflection 
field will change and therefore also the surface density and 
the magnification factors. In particular, the relation in equa- 
tion (|20|l will no longer hold. However, since this no longer 
corresponds to a rescaling of the sources, it can hardly be 
considered as caused by the mass sheet degeneracy. 



5 EXPERIMENT 

5.1 Positional information 

To further illustrate the significant practical problems posed 
by these degeneracies, we shall use a simulated lensing sce- 
nario. The left panel of Fig.|6]shows the true mass map used 
in the simula tion, an elliptical generalization of a NFW mas s 
distribution (|Navarro et al. I ll996l ; IWright fc Brainerdll2000T l , 
with concentration parameter c v i r = 5 and scale parameter 
r s = 50 arcseconds. The lens was placed at a redshift of 0.7 
in a cosmology with Ho = 70 km s -1 Mpc -1 , f2 m = 0.3 
and S1a = 0.7. The surface density is expressed in terms of 
the critical density for a source at z = 3. Eight sources were 
placed at different redshifts. These generate 26 images which 
can be seen in the contour map, and also in the right panel 
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Figure 6. Left panel: the input mass distribution for a simulated lensing scenario consists of an elliptical generalization of a NFW profile. 
Eight distinct sources generate 26 images which are shown on top of the contour map. Center panel: for a source at z = 3, the mass 
distribution from the left panel generates these critical lines and caustics. Right panel: the same 26 input images as in the left panel are 
shown here. The three images of a single source which will be used for additional time delay constraints later, are surrounded by squares. 




Figure 7. Left panel: when only positional inform ation about the ima ge systems is used and only Plummer basis functions are used, 
corresponding to the algorithm as described in lLiesenborgs et al. the average of 20 individual solutions yields this mass distribu- 

tion. Although the overall shape is similar to the mass distribution in Fig. [6] this mass map is clearly much more steep. Center panel: the 
critical lines and caustics for this average solution. Comparing this to the center panel of Fig. [6] one immediately notices that the scale of 
the caustics is much larger, indicating the presence of a mass sheet-like degeneracy. Right panel: the densities of the individual solutions 
are compared to the true input mass densities, precisely at the locations of the image centers. Being much more steep and apparently 
lacking a constant offset again reveals the presence of a mass sheet-like degeneracy. 



of the same figure. The redshifts of the sources lie between 
z = 2.7 and z — 3.4. The center panel of Fig. [6] shows the 
critical lines and caustics for a source at z = 3. We shall 
be interested in comparing the scale of the caustics here to 
that of reconstructions, as rescaling the source plane is the 
hallmark of the mass sheet degeneracy. 

First, a reconstruction i s done using the algorithm pre- 
sente d in Liesenborgs et alj |2006). In this method, Plum- 
mer (|Plummerlll91ll ') basis functions are arranged according 
to an adaptive grid, and the weights of these basis functions 
are determined using a genetic algorithm. The left panel of 
Fig. [JJ shows the resulting mass distribution, an average of 
20 individual reconstructions. Overall, it has a similar ap- 
pearance as the input mass map shown in Fig. [(J and its 
critical lines (center panel) do bear resemblance to those of 
the true lens. The resulting mass distribution, however, is 
considerably steeper than the one used in the simulation, 
which can also clearly be seen in the right panel of Fig. [JJ 
In this image we show the densities measured at the cen- 
ter of each image. For the sake of visualization, instead of 



using the two dimensional coordinate of each image center 
point, the distance to the NFW center is used on the a;-axis. 
Comparing the two profiles on that figure, one is immedi- 
ately reminded of the mass sheet degeneracy, i.e. the profiles 
seem to differ mostly by a rescaling and the addition of a 
specific offset. Similarly, when looking at the caustics in the 
center panel of Fig. [JJ they seem to be an expanded version 
of the ones of the input lens, again a feature that hints at 
the presence of the mass sheet degeneracy. 

Due to the inherent uncertainty in between the images 
caused by the monopole degeneracy, the most sensible thing 
to do in trying to retrieve the NFW parameters, is to per- 
form a fit based on the densities precisely at the locations of 
the images. When looking at the steepness of the recon- 
structed mass map, it will come as no surprise that ex- 
tracting the concentration and scale parameters from the 
reconstructed densities will not yield accurate results. In 
this particular case, the estimates are cvir = 14.1 jig's and 
r B — 11.71q'7 arcseconds respectively (68% confidence level). 

Using Plummer basis functions it is of course rather dif- 
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Figure 8. Loft panel: to generate this aver age solution, the same alg orithm was used as in Fig. [7] but this time a constant density sheet 
was allowed as a basis function as well (sec Liesenborgs ct al . (2009)). This solution resembles the input mass distribution more closely. 
Center panel: when the critical lines and caustics for this solution are calculated, the source plane no longer exhibits an obvious rescaling 
compared to the input situation. Right panel: comparing the densities at the image centers to those of the input mass map confirms that 
this reconstruction is much more similar to the true situation than the mass distribution from Fig. [7] 



ficult to mimic the effect of a sheet of mass. Given the fact 
that within the images a considerable mass sheet-like com- 
ponent seems to be present in the input lens, the fact that 
the reconstruction is rather different is not that surprising. 
In a second reconstruction not only Plummer basis functions 
were used, but a single mass s heet basis function was p resent 
as well, as is also explained in lLiesenborgs et alJ (|2009h . The 
resulting mass distribution, again an average of 20 individual 
solutions, can be seen in the left panel of Fig. [8] This mass 
distribution is clearly a lot more similar to the input mass 
map, also obvious from the right panel of the same figure. 
As can be expected from the similarity between the true and 
reconstructed lenses, the caustics shown in the center panel 
no longer show a clear difference in scale. When performing 
a fit in a similar fashion as the previous reconstruction, one 
finds NFW parameters of c vir = 5.9±o;g and r s = 38.3±4'.4 
arcseconds, corresponding better to the true values. 

Unfortunately, in this example there is no real reason 
to prefer one reconstruction over the other, meaning that 
apart from the overall shape of the true mass distribution, 
little would be revealed by performing these inversions. Fur- 
thermore, in both cases fitting the NFW shape yields an 
overestimate of the concentration parameter and an under- 
estimate of the scale radius, although the effect is far more 
serious in the first inversion. 



5.2 Including time delay information 

Earlier in this article it was mentioned that time delay in- 
formation can help not only in breaking the most basic mass 
sheet degeneracy; also for the generalized versions time de- 
lays depend on the precise construction meth od of the degen- 
erate mass map. In lLiesenborgs et al. | l|200gh a method was 
described with which time delay information can be added 
as a constraint, now using a so-called multi-objective genetic 
algorithm. To explore the value of such time delays, time 
delay information of a single system of images was added 
to the inversion procedure. The three images used for this 
added constraint are indicated in the right panel of Fig. [5] 
and correspond to a lensed source at redshift 2.75. 

First, the inversion was run using Plummer basis func- 



tions alone, i.e. no mass sheet basis function was used. The 
resulting mass distribution can be seen in the left panel of 
Fig. [9] Concerning the input to the inversion procedure, the 
only difference with the results obtained in Fig.[7]is the addi- 
tion of time delay information. The effect of this information 
is impressive however: especially when comparing the right 
panels of these two figures it is clear that the time delay 
constraint seems to pull the densities at all image locations 
towards the true densities. The caustics in the center panel 
confirm the mass sheet-breaking effect: the overall scale of 
the source plane seems to match the true one more closely. 
Note that this is solely the effect of time delays for a single 
image system. When determining the NFW parameters as 
before, one now finds c v i r = g and r B = 28.8 j^f arcsec- 

onds. This result still does not reflect the true parameters, 
but is far less deviant than the estimates based on Fig. (7) 

As a final inversion, both the time delay information 
and the mass sheet basis function was used, i.e. the inversion 
is similar to the one reflected by Fig. [8] but now including 
time delay information about one set of images. The result- 
ing mass map is depicted in the left panel of Fig. 1101 the den- 
sities at the image locations can be seen in the right panel 
of the same figure. The effect of the added time delay infor- 
mation is not as impressive as in the previous example, but 
still it leads to an improvement of the estimated densities at 
most image points. Only the density at the point closest to 
the NFW center seems to be underestimated systematically. 
This single outlier does not seem to have a profound adverse 
effect on the estimation of NFW parameters however, which 
now are c v - 1T = 4.8^0* and r s — 52.1^4'j arcseconds, very 
close to the true values. 



6 DISCUSSION AND CONCLUSION 

The examples shown in this article make it clear that these 
lensing degeneracies are closely related to substructure in 
the mass density. The monopole degeneracy is the worst, 
as it has no effect on any property of the images: the same 
sources predict the same images, with identical magnifica- 
tions and identical time delays. In essence it means that we 
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Figure 9. Left panel: this is the average solution when only Plummer basis functions are used, similar to Fig. [3 but as an additional 
constraint the time delays between the images indicated in the right panel of Fig. [6] were included. Comparing this result to the one in 
Fig. [7] reveals that the addition of this information has a profound effect on the steepness of the solution. Center panel: thanks to the 
presence of the time delay information the caustics of this average solution no longer display the obvious rescaling that could be observed 
earlier. Right panel: as was already suggested by the left panel, comparing the densities at the image locations shows that the time delay 
information for even a single image system seems to pull all reconstructed densities towards the true densities. 
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Figure 10. Left panel: when both time delay information about the system marked in the right panel of Fig. [6] and a mass sheet basis 
function are allowed, this average result is obtained. Apart from the sharp central peak, the mass distribution is very similar to the input 
mass map. Center panel: calculating the critical lines and caustics of this solution reveals that the source plane is very similar in scale 
to the one of the input lens. Right panel: in this case the densities at the image locations are extremely similar to the input densities. 
Only the density at the point closest to the center of the input NFW distribution appears to be underestimated systematically. 



do not have a firm handle on the mass distribution in be- 
tween the images, the only real constraint being the absence 
of unobserved images. This illustrates both the importance 
of gravitational lens systems with a multitude of images and 
when necessary, of the use of the null space, i.e. the locations 
where no images are observed. 

In its most general version, the mass sheet degeneracy 
does not require an actual sheet of mass, nor does it neces- 
sarily modify the slope of the mass distribution, as the al- 
ternative name of steepness degeneracy might suggest. Com- 
paring equations (|17[) and (|20[1 . it becomes clear that 
when scaling source A with factor Aa and source B with 
factor Ab, the effect on the mass density is the following: 

Ei (J A ) = A a So (/a) + (1 - A A )E cr (z A ) 

Ex(/b) = A B Eo (J b ) + (1 - A B )E cr (zb) , (21) 

where 7a and Ib are again the locations of the images of 
source A and B respectively. Of course, this can be gen- 
eralized to any number of sources and images. Apart from 
rescaling each source in the system, equation (|21|l illustrates 



that this degeneracy changes the density precisely at the 
locations of the images, i.e. it necessarily introduces sub- 
structure. The effect on the time delays is in general not 
straightforward implying that time delay information may 
be very valuable in helping to break this degeneracy. The 
simple effect on the magnification is not easy to use since 
one most often does not have information about the size 
of the source. Efforts to use the magnification information 
to avoid the m ass sheet deg e nerac y have been undertaken 
however, e.g. in lTavlor et al.l () 1998T) . 

The fact that both monopole and mass sheet degenera- 
cies are - at least in principle - possible with any amount 
of sources and images, combined with the fact that they are 
both linked to substructure in the mass distribution, means 
that in reality it is the prior information on the degree of 
smoothness of the mass distribution that effectively breaks 
these degeneracies. The degree of smoothness can depend 
implicitly on the method used, e.g. in para metric methods 
like the Len s Tool method (I Julio et al.ll2007t ), the method of 
IZitrin et all (|2010h which is largely based on the distribution 
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of the light, or the method of the authors in which overlap- 
ping s mooth basis functions are used (e.g. iLiesenborgs et all 
(2007)). Other inversion procedures use an explicit regular- 
ization scheme or an explicit Bayesian prior, e. g . in th e Pix- 
eLens method (|Saha fc Williams! l|2004 ), IColesl (|2008l V). The 
dependence on explicit or implicit assumptions about the 
smoothness of the solution of course means that one must 
be careful when interpreting obtained inversion results, not 
only in lensing systems with few sources, but also when a 
larger amount of sources are available, as indicated by the 
experiment. This same experiment also highlights the im- 
portance of time delay information, which easily asserts a 
global effect on the reconstruction. It is unfortunate that in 
many cases of interest these time delays may simply be too 
long to measure in practice. 



This paper has been typeset from a Tp^X/ DTffjX file prepared 
by the author. 
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